Theory of laser-induced demagnetization at high temperatures 



A. Manchon^'^^ lji^ Xu\ and S. Zhang^ 
^Department of Physics, University of Arizona, Tucson, AZ 85721, USA; 
"^Materials Science and Engineering, Physical Science and Engineering Division, KAUST, Saudi Arabia 

(Dated: December 13, 2011) 

Laser-induced demagnetization is theoretically studied by explicitly taking into account interac- 
tions among electrons, spins and lattice. Assuming that the demagnetization processes take place 
during the thermalization of the sub-systems, the temperature dynamics is given by the energy 
transfer between the thermalized interacting baths. These energy transfers are accounted for explic- 
itly through electron-magnons and electron-phonons interaction, which govern the demagnetization 
time scale. By properly treating the spin system in a self-consistent random phase approximation, 
we derive magnetization dynamic equations for a broad range of temperature. The dependence of 
demagnetization on the temperature and pumping laser intensity is calculated in detail. In par- 
ticular, we show several salient features for understanding magnetization dynamics near the Curie 
temperature. While the critical slowdown in dynamics occurs, we find that an external magnetic 
field can restore the fast dynamics. We discuss the implication of the fast dynamics in the application 
of heat assisted magnetic recording. 
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I. INTRODUCTION 

Laser-induced Dcmagnctizationii— (LID) and Heat As- 
sisted Magnetization Reversal^ (HAMR) constitute a 
promising way to manipulate the magnetization direc- 
tion by optical means. While both LID and HAMR in- 
volve laser-induced magnetization dynamics of magnetic 
materials, there are several important differences. LID 
is usually considered as an ultrafast process vifhere the 
hot electrons excited by the laser field transfer their en- 
ergy to the spin system, causing demagnization. The 
demagnetization time scale ranges from 100 femtosecond 
to a few picoseconds. For HAMR, the laser field is to 
heat the magnetic material up to the Curie temperature 
so that the large room-temperature magnetic anisotropy 
is reduced to a much smaller value and consequently, a 
moderate magnetic field is able to reverse the magneti- 
zation. The time scale for the HAMR process is about 
sub-nanosecond, three orders of magnitude larger com- 
pared to LID. 

LID observations have been carried out in a 
number of magnetic materials including transition 
metalsii^"— , insulators^, half-metals^rJ^ and dilute mag- 
netic semiconductorsii. A general consensus of the laser- 
induced demagnetization process is that the high energy 
non-thermal electrons generated by a laser field relax 
their energy to various low excitation states of the elec- 
tron, spin and lattice^^. The phenomenological model for 
this physical picture is referred to as three-temperature 
model^i^ where the three interacting sub-systems (elec- 
trons, spins, lattice) are assumed thermalized individ- 
ually at different temperatures which are equilibrated 
according to a set of energy rate equations. By fitting 
experimental data to the model, reasonable relaxation 
times of the order of several hundred femtosecond to a 
few picoseconds have been determined. 



Various microscopic theories^'^^^— have been proposed 
to interpret these ultrafast time scales of electron-spin 
and electron-lattice relaxations. Zhang and Hiibner— 
proposed that the laser field can directly excite the spin- 
polarized ground states to spin-unpolarized excited states 
in the presence of spin-orbit coupling, i.e., the spin-flip 
transition leads to the demagnetization during the laser 
pulse. In this picture, the demagnetization is instan- 
taneous («50-150fs). Recent numerical simulationsii 
show that due to a few active "hot spots", the instan- 
taneous demagnetization is expected for at most a few 
percent of the magnetization, consistently with experi- 
mental argumentsi^. Koopmans et ali^i^ suggested that 
the excited electrons lose their spins in the presence of 
spin-orbit coupling and impurities or phonons, through 
an " Elliot- Yafet" -type (EY) spin- flip scattering. Recent 
numerical evaluations of the EY mechanism in transition 
metalsi^ tend to support this point of view. Alterna- 
tively, Battiato et alJ^ recently modelled such ultrafast 
demagnetization in terms of superdiffusive currents. Fi- 
nally, numerical simulations of the ultrafast demagneti- 
zation based on the phenomenological Landau-Lifshitz- 
Bloch equation have been achieved successfully^. 

While these demagnetization mechanisms provide rea- 
sonable estimation for the demagnetization time scales, 
the theories are usually limited to the temperature much 
lower than the Curie temperature and/or make no di- 
rect connection to the highly successful phenomenologi- 
cal three-temperature modcli~— . As it has been recently 
shown experimentallyii^, most interesting magnetization 
dynamics occur near the Curie temperature. 

In this paper, we propose a microscopic theory of the 
laser-induced magnetization dynamics under the three- 
temperature framework and derive the equations that 
govern the demagnetization at arbitrary temperatures. 
More specifically we predict magnetization dynamics in 
the critical region. 
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The paper is organized as follows. In Sec. |TT1 we pro- 
pose a model for LID processes. In Sec. Illli we describe 
the spin system by the Heisenberg model which is solved 
by using a self-consistent random phase approximation. 
In Sec. lIVi the central dynamic equations for the magne- 
tization are derived. In Sec.|Vl the numerical solutions of 
the equations are carried out and the connection of our 
results with the experimental data of LID and HAMR is 
made in Sec. I VII We conclude our paper in Sec. I VIII 

II. MODEL OF LID 
A. Spin loss mechanisms 

One of the keys to understand ultrafast demagneti- 
zation is to identify the mechanisms responsible for the 
spin memory loss. In the case of transition metal ferro- 
magnets for example, the spin relaxation processes lead 
to complex spin dynamics due to the itinerant character 
of the magnetization. Elliott^ first proposed that dclo- 
calized electrons in spin-orbit coupled bands may lose 
their spin under spin-independent momentum scatter- 
ing events (such as electron-electron or electron-impurity 
interaction). This mechanism was later extended to 
electron-phonon scattering by Yafet and Overhauser— . 
Consequently, the spin relaxation time Tg is directly pro- 
portional to the momentum relaxation time Tp. Whereas 
the electron-electron relaxation time is on the order of 
a few femtoseconds^^ (fs), the electron-impurity and 
electron-phonon relaxation time is on the picoseconds 
(ps) scale. In semiconductors, bulk and structural inver- 
sion symmetry breaking as well as electron-hole interac- 
tions lead to supplementary spin relaxation mechanisms 
such as D'yakonov-Perel2i and Bir-Aronov-Pikus^^ that 
are beyond the scope of the present study. 

Relaxation processes also apply to collective spin exci- 
tations such as magnons. Whereas the electron-magnon 
interaction conserves the angular momentum, magnon- 
magnon interactions and magnon-lattice interactions in 
the presence of spin-orbit coupling contribute to the total 
spin relaxation. While the former occurs on the magnon 
thermalization time scale^^ (lOOfs), the latter is however 
at the second order in spin-orbit coupling and is con- 
sidered to occur on the lOOps time scale. Therefore, in 
a laser- induced demagnetization experiment, it is most 
probable that all the processes mentioned above take 
place during the thermalization time scale of the excited 
electrons and excited magnons. 

B. Demagnetization scenario 

To establish our model, we first separate the LID pro- 
cesses into four steps: (i) generation of non-thermal hot 
electrons by laser pumping; (ii) relaxation of these hot 
electrons into thcrmalized electrons characterized by an 
electron temperature Tg; (iii) energy transfer from the 



thermalized hot electrons to the spin and lattice sub- 
systems; (iv) heat diffusion to the environment. 

In our model, to be given below, we will take steps 
(i) and (ii) infinitely fast. In the step (i), a laser pump 
excites a fraction of electrons below the Fermi sea to ssl.5 
cV above the Fermi level. This excitation process is of the 
order of a few fs. The photo-induced electron transition is 
considered spin conserving and thus does not significantly 
contribute to the demagnetization although the spin-flip 
electron transition could occur in the presence of the spin- 
orbit couplingi^. 

In step (ii), the strong Coulomb interaction among 
electrons relaxes these non-thermal high-energy electrons 
to form a hot electron bath which may be described by 
a thermalized hot electron temperature Tg. During this 
electron thermalization process, strong electron-electron 
interaction-induced momentum scattering in the pres- 
ence of spin-orbit coupling leads to the ultrafast transfer 
of the spin degree of freedom to the orbital onoSi. In our 
model, the electron thermalization is considered instan- 
taneous and any possible femtosecond coherent processes 
are disregarded^^. Therefore, due to ultrafast (fs) mo- 
mentum scattering, the thermalized hot electrons act as 
a spin sink. Under this approximation, the demagnetiza- 
tion itself, defined as the loss of spin angular momentum, 
takes place during the thermalization of the electron bath 
in the presence of (either intrinsic or extrinsic) spin-orbit 
coupling. 

Following the definition of the three-temperature 
model, we assume that the system can be described 
in term of three interacting baths composed of laser- 
induced hot electrons, spin excitations of the ground 
state (magnons) and lattice excitations (phonons). The 
applicability of this assumption is discussed in Sec. IIIDI 
Therefore, the magnetic signal essentially comes from the 
collective spin excitation and it is assumed that the laser- 
induced hot electron only contribute weakly to the mag- 
netization. Consequently, under the assumption that the 
spin loss occurs during the thermalization time of the 
electron and spin systems, the demagnetization problem 
reduces to tracking the energy transfer between the spin 
bath and the electron and phonon bathd. 

Our main objective is then to understand step (iii), 
where the electrons at a higher temperature transfer 
their energy to the spin and lattice sub-systems. Un- 
der the electron-magnon interaction, the magnons spin is 
transferred to the electron system, and is eventually lost 
through thermalization of the electron bath. Through 
interactions among electrons, spins and lattice, the en- 
tire system will ultimately reach a common temperature. 
Finally, a heat diffusion, step (iv), will expel the heat to 
the environment; this last step will be considered via a 
simple phcnomcnological heat diffusion equation. 

To quantitatively determine the energy transfer among 
electrons, spins and lattice in the step (iii), one not only 
needs to know the explicit interaction, but also the distri- 
bution of the densities of excitations (electrons, magnons 
and phonons). Within the spirit of the three tempera- 
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ture model, we consider that each sub-system (electron, 
spin and lattice) is thermalized, i.e., one can define three 
temperatures for electrons Te, spins Tg and lattice T;. 
The justification of this important assumption has been 
made in the previous section and can be qualitatively 
summarized: 1) For the hot electrons of the order of 
leV, the electron-electron relaxation time is T^e ~ 10/s, 
which is about 100 times faster than the electron-spin and 
electron-phonon interactions^^. 2) The lattice-lattice in- 
teraction is about one order of magnitude smaller than 
the electron-electron relaxation time, tu w lOO/s^S.. 3) 
Multiple spin-waves processes are known to take place 
in the ferromagnetic relaxation leading to so-called Suhl 
instabilities^. The relaxation time is of the order of 
Tss oc Ti/Tc ~ 100/s at least for high energy magnons^^ 
(for long wave length magnons, the lifetime could be 
significantly longer). Thus, it is reasonable to assume 
that the concepts of the three temperatures are approx- 
imately valid as long as the time scale is longer than 
sub-picoseconds. 



C. Model Hamiltonian 

We now propose the following Hamiltonian for LID 

H = J2h^ + Hes + H,i + Hsi : (1) 



where {fi = e, s, /) are the electron, spin and lat- 
tice Hamiltonians, and H^i, {fj, ^ v) are the interaction 
among sub-systems. In the remaining of the article, the 
hat ' denotes an operator. Each term is explicitly de- 
scribed below. 

The electron system is described by a free electron 
model He = X^k ^kcj^Ck where cj^ (ck) represents the 
electron creation (annihilation) operator. The equilib- 
rium distribution is simply the Fermi distribution at Tg. 
The lattice Hamiltonian Hi = J2qx ^'^kA^kA^kA is mod- 
eled by simple harmonic oscillators where b^^ {bk\) is the 
phonon creation (annihilation) operator and A is the po- 
larization of the phonon. The phonon distribution at T/ 
is «kA = [exp(7ia;^^/fcs21) — The spin Hamiltonian 

is modeled by the Heisenberg exchange interaction. 



H. 



— — ^ JijSi ■ Sj — g^sHe 



(2) 



<i'j> 



where Jy is the symmetric exchange integral. Si is the 
spin operator at the site and H^x is the external mag- 
netic field applied in z-direction. Unlike the electron and 
lattice Hamiltonians, the spin Hamiltonian is not a single 
particle Hamiltonian and the distribution of the spin den- 
sity is neither a fermionic nor a bosonic distribution. To 
describe the equilibrium distribution of the spin system 
at arbitrary temperatures, we will model the equilibrium 
properties of the spin system in the next section. 



The electron-lattice interaction H^i is taken as a stan- 
dard formes 

Hei = ^Bqx{ci+^cJ}qx + c+_^cJ)+^), (3) 

k,qA 

where the Bc^x is the electron-phonon coupling constant. 
For acoustic phonons, the coupling constant takes a par- 
ticularly simple formes. 



SqA — 



^Y2MA^- 



(4) 



Here cf is the electron Fermi energy and M is the mass 
of the ion. 

The electron-spin interaction H^s is modeled by the 
conventional exchange interaction (sd Hamiltonian) : 



H, 



-Jex c^e^"" ""' ■ SjOck'e-*''' "-^ , (5) 

i,k,k' 



where we have assumed a constant coupling constant Jex 
and (T is the electron spin. When one replaces a ■ Sj by 
&zSj + i(tT_S'^ -I- a+S~), the above Hes contains two 
effects: the first term is responsible for the spin-splitting 
of the conduction bands and the second term leads to a 
transfer of angular momentum between the spins of the 
hot electrons and the spins of the ground state, i.e. spin- 
waves generation and annihilation. While the interaction 
conserves the total spin angular momentum, the thcr- 
malization process of each bath is not spin conserving 
as mentioned above. Therefore, this interaction trans- 
fers energy between the electron and spin baths, which 
results in the effective demagnetization of the magnon 
bath. Consequently, the generation of magnons by hot 
electron is a key mechanism in our model (see also Rcf. 

Finally, the spin-lattice interaction Hsi has been at- 
tributed to spin-orbit coupling^. The energy and the 
angular momentum conservations require Hsi containing 
two-magnon (flqCiq') and two-phonon operators (&j^^k')- 
Since the spin-orbit coupling is already treated as a 
perturbation, this process is second order in the spin- 
orbit coupling parameter and it is expected to be rather 
small^. Thus, Hsi is much smaller than Hes and Hei, 
and we place Hsi = throughout the rest of the paper. 

To summarize our model, we consider three subsystems 
(electrons, spins, and lattice) described by He, Hs and 
Hi respectively. These subsystems have their individual 
equilibrium temperatures Te, Ts and T/. The heat or 
energy transfer among them are given by the interaction 
Hes and Hei ■ To determine the kinetic equation for three 
subsystems, we should first establish the low excitation 
properties of the spin system from Hg and relate Ts to 
the magnetization 'm[Ts). 
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D. Materials considerations 

As stated in the introduction, laser-induced demagne- 
tization has been observed in a wide variety of materials 
presenting very diverse band structures and magnetism. 
From the materials viewpoint, the present model makes 
three important assumptions: (i) laser-induced hot elec- 
trons, ground state spin excitations and phonons can be 
treated as separate interacting sub-systems; (ii) there ex- 
ists a direct interaction between hot electrons and collec- 
tive spin excitations; (iii) the excited spin sub-system can 
be described in terms of spin-waves. 

Whereas the consideration of a separate phonon bath 
is common, the separation between the electron and spin 
populations may seem questionable. In systems where 
the itinerant and localized electrons can be identified 
(such as 4/-rare earth or carrier-mediated dilute mag- 
netic semiconductors), it seems quite reasonable. How- 
ever, in typical itinerant ferromagnets such as transi- 
tion metals, the magnetism arises from a significant por- 
tion of itinerant electrons. We stress out that in our 
model, the separation between electron and spin baths 
arises from the fact the electrons we consider are laser- 
induced hot electrons near Fermi level (in the range 
[ep — ksTe, ep + fcs?e]), whereas the spin bath describes 
the magnetic behavior of electrons lying well below Fermi 
level. The concept of spin waves used in the present ar- 
ticle is rather general and applies to a wide range of fer- 
romagnetic materials. Although energy dispersion may 
vary from one material to another, it is unlikely to have 
strong influence on the main conclusions of this work. 

The interaction between hot electrons and magnons is 
actually more restrictive since it assumes overlap between 
electrons near and far below Fermi level. For example, 
this approach does not apply to half-metals (electron- 
magnon interaction is quenched by the 100% spin polar- 
ization) or magnetic insulators. Nevertheless, in metallic 
materials such as transition metals and rare-earth, this 
interaction does not vanish and can lead to strong spin 
wave generation, as demonstrated by Schmidt et al*2i in 
Fe. 



III. EQUILIBRIUM PROPERTIES OF THE 
SPIN SYSTEM 



scribe below a self-consistent random phase approxima- 
tion which reproduces Bloch's law at low temperatures 
and the mean field result at high temperatures. 

We first recall some elementary relations of these spin 
operators given below. 



qX ■ qy 



2Sf 



(6) 
(7) 

S+S- ^S{S + l) + S!-{St)\ (8) 
and the spin Hamiltonian, Eq. ^ can be rewritten as 

ij i 

Our self-consistent random phase approximation treats 



the resulting commutator. 



2Sf 



2m{T) 



c-number, where m(T) is the thermal average of 
Sf to be determined self-consistently. If we introduce 
the Fourier transformation, ^ (1/A^) •^i'^e"*-^', 
the above commutator reads as [5*^,5*"] = 2m(r)Jkq 



and thus by introducing dj^ 
imuti 

Similarly, we have [Hg , d, 



Sf/y/2m{T), one has 



a standard boson commutator relation [dq, d^, 

TiWqdq, where 



huq^ gHBHecc + 2m{T)J2[Jo-Jiq)] (10) 

q 

where J(q) = {l/N)J2<ij> Ji] cxp[iq ■ (R^ - R^)]. 
With the above bosonic approximation, one can self- 
consistently determine the magnetization m(T) and 
other macroscopic variables such as the spin energy and 
specific heat. A particular simple case is for the spin-half 
= 1/2 where the identity 



S, = S- S-S+ = 1/2 - E 2m(T)a+( 



(11) 



immediately leads to the self-consistent determination for 
m(T) 



(1^^ 1/9 1 2m(r) 
m{T) = 1/2 ---^-^j—^ 



(12) 



The Heisenberg model for the spin system, Eq. ([2|), 
has no exact solution even in equilibrium. At low tem- 
perature, the simplest approach is based on the spin- 
wave approximation which predicts Bloch's law for the 
magnetization m(T) = mp — B{T/Tc)^/^ where Tc is 
the Curie temperature and B is a numerical constant^^. 
As the temperature approaches the Curie temperature, 
Bloch's law fails. Instead, one uses a molecular mean 
field to model the magnetization. The resulting magne- 
tization displays a critical relation near Tc, i.e., m(T) oc 
(1 — T/TcY/"^. Since we are interested in modeling the 
magnetization in the entire range of temperature, we de- 



At low temperature, one can approximately replace Tn{T) 
by 1/2 in the right-hand side of the equation and one 
immediately sees that the above solution produces the 
well-known Bloch relation, i.e., 1/2 — m(r) cx T^/^. Near 
the Curie temperature, one expands e^^'^'^ = 1-1- /3?iajq 4- 
(l/2)(/3?ia;q)^ and notice that cjq is proportional to ra{T) 
at zero magnetic field, see Eq. pUj) . By placing this 
expansion into Eq. the zero order term in m(T) 

determines the Curie temperature and the second order 
term gives the scaling rri^iT) oc {Tc — T), i.e. the mean 
field rcsuh is recovered, m{T) cx (1 - T/Tc)^/'^. Thus 
the self-consistent approach captures both low and high 
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temperature limiting cases. In fact, the Green's function 
technique^^ has been developed to justify this approxi- 
mation. 

For the cases other than S = 1/2, the relation between 
Sf and the number of magnons is more complicated due 
to non-constant (5f )^ and thus Eq. ([8]) cannot immedi- 
ately lead to a self-consistent equation for m(T). Instead, 
one needs to relate ((^f )^) to m(T) and the magnon den- 
sity. TyablikoA*^ introduces a decoupling method to ap- 
proximate ((S'f )^) with ni{T) and the normalized number 
of magnons 



no 



exp(/?wq) - 1 



(13) 



Here finds that, for arbitraryS", the self-consistent equa- 
tion for determining Tn(T) is 

_ {S- no){l + n,fS+\l + S + no)nf +^ 



,2S+1 
*0 



By replacing S — 1/2. the above equation reduces 
to Eq. (jl2p . The magnetic energy can be similarly 



2ni 







exp(/3a;q) - 1 



where Eq is the ground state energy and hujg{0) is the 
spin wave energy at T = 0. Once the internal energy is 
obtained, the specific heat, Cp ~ dE/dT, may be numer- 
ically calculated. 

m{Ts) is uniquely determined from Eq. (14) or Eq. (12) 
for s = 1/2, if the spin temperature is known. Thus, the 
laser-induced demagnetization is solely dependent on the 
the time-dependent spin temperature Tg . Before we pro- 
ceed to calculate Ts{t) or m{t), we show the solutions of 
Eq. (14) or Eq. (12). In Figure [U the reduced magneti- 
zation m{T)/ S and the specific heat as a function of the 
normalized temperature T/T^ with [Figs. [TJa) and (b)] 
and without [Figs. [iKc) and (d)] the magnetic field are 
shown. A few general features can be readily identified. 
First, the shapes of the magnetization curves for differ- 
ent spins are very similar. Second, the magnetic field 
removes the divergence of the specific heat at the Curie 
temperature. As expected, the magnetization reduces to 
that of the mean field result near the Curie temperature 
and to that of the spin wave approximation at low tem- 
peratures. 



IV. DYNAMIC EQUATIONS 

The energy or heat transfer among electrons, spins and 
lattice may be captured by the general rate equations 
given below, 



d£e 
dt 



(16) 




FIG. 1: (Color online) Temperature dependence of (a) 
magnetization and (b) specific heat (arbitrary unit) for 
spin=l/2,l,2,8 in the absence of the external field; Temper- 
ature dependence of (c) magnetization and (d) specific heat 
for spin=l/2,l,2,8 in an external field H/Tc = 0.001. 



-TT — ^ si + -L eZ 

at 

—1 = r 

dt 



(17) 
(18) 



where £i are the energy densities {i = e, s,l) and the 
rate of the energy transfer Tij should be determined by 
Eq. (1). Since we have neglected the weaker interaction 
between spins and lattice, we set T^i = in the above 
equations. In the following, we explicitly derive the re- 
laxation rates of Fe; and T^s from Eqs. (3) and (5). 



A. Electron-lattice relaxation T^i 

The energy transfer rate between electrons and 
phonons does not involve the spin. The Fermi golden 
rule applied to Eq. (3) immediately leads to 



47r 



(19) 



k,q 



(71k+q(l - ?ik)(l + J^q) - ?^k(l - nk+q)nP), 

where the first (second) term represents the energy trans- 
fer from (to) the electrons to (from) lattice by emit- 
ting (absorbing) a phonon. Note that the electrons and 
phonons have different temperatures; otherwise the de- 
tailed balance will make the net energy transfer zero. The 
electron and phonon densities are given by their respec- 
tive equilibrium temperatures at Tg and Tj, i.e., = 
[exp((ek-ef)/fcBTe) + l]-i and = [exp(nu;P//cBTi) - 
1]~^. We consider polarization-independent acoustic 



phonons, i.e. 



VsQ where Vs is the phonon veloc- 



ity. By replacing Bp given in Eq. (4) into Eq. (19), we 
have 
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47r /2 y V merrie 3 - /■+°° e^^de 

where we have defined the cut-off energy Cq = (q — which comes from the 5 function, and we have 

introduced the maximum phonon wave number in the First Brillouin zone (this definition is the same for magnons 
and Fermi wave vectors), ~ kp — (67r^)^/"^/ao. By integrating over the electron energy e, we obtain 

Ait [2 V me nie f"" 3,_v_„o/rT.N .^v/rT.^^-|^ rj. I ^^'sQ 



Tei^-i^ep] —-^^ / g^dgK(re) - n^(TO)fcBTe ^ - In I t:^^^ | | , (21) 




where we have defined a Debye temperature, 9 = 
JiVsQm/kB- The last term of Eq. (j2ip can be approxi- 
mated by 




hvsq 



Q{2kF~p), 
(22) 



for hujq <^ ksTe, where 8(x) is the step function. There- 
fore, the relaxation rate becomes 



1 

a V 3 



Ga 



Te 



Ga 



Ti 



and Gn{,x) 



(23) 



/;'^fW(e*-i). 



Interestingly, for Te,Ti > 0, the relaxation rate reduces 



to: 



1/2 \ 97r me /T,. - T/ 



(24) 



Thus, the relaxation rate is simply proportional 
to the difference between the electron and lattice 
temperatures (Fej (x Tg — Tj); this is the assumption made 
in the earlier three temperature model^. 



B. Electron-spin relaxation Fes 

The interaction between the electrons and spins given 
by Eq. (5) may be simplified by using the self-consistent 
random phase approximation, i.e., we replace S'f by its 
thermal average m{T,) and = y'2m{Ts)a^. The 
electron-spin interaction can then be rewritten as 



J e: 



kq 



(25) 



where we have dropped the m{T)az term since it does 
not involve the energy transfer between the electron and 
the spin. 



The second order perturbation immediately leads to 
the electron-spin relaxation: 

2TT2Sm{Ts) j2 \r% p. \ 



TV 



k,q 



X (nk;(l - nk-qt)(l + Wq) - nk-qt(l - "k;)'^q) 

(26) 

where Wq is the magnon frequency given by Eq. (10), the 
electron distribution is riko- = [cxp((ek — ei?)/fcBre) + l]~^ 
and the magnon distribution is riq — [exp(?ia;q) — 1]^^. 
Note that the electron sub-system is considered unpolar- 
ized due to the strong spin relaxation occuring during 
thermalization. In Eq. (26), the first (second) term rep- 
resents the electron emitting (absorbing) a magnon. For 
the long wavelength, the magnon dispersion, Eq. (10), is 
simply hujq — pLsHext + otkBTcq^a^ where awl. Fol- 
lowing the same procedure as the previous section, we 
find, 



— 2SM{T)J^ ^ — 



X / qdq{nuJqY{n'q{Te)-n'^{T,)) 
10 



(27) 



If the magnetic field is zero, the integration over q can 
be immediately carried out and we approximately have 
(discard the numerical constant a) 



f Ts 



2hV 

DTr 



G2 



(28) 



where Gn{x) has been defined below Eq. (23) and the 
temperature-dependent spin stiffness is _D = ra{Ts)q^a? . 
An important result of this paper is that Fe^ is propor- 
tional to rrv^iTs) and thus it is vanishingly small near 
the Curie temperature. Furthermore, since Tg/Tc and 
Ts/Tc are always comparable to 1, the electron-spin re- 
laxation rate given by Eq. (28) is not proportional to 
Te — Ts, which is quite different from the previous three- 
temperature modeli. 
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C. Specific lieats of subsystems 

As the right sides of the rate equations, Eqs. (16)-(18), 
are expressed in terms of three temperatures Tg, and 
T/, we need to relate the energy change of each system 
to the temperatures, i.e., we should define the heat ca- 
pacities Ci for each subsystem dEi = CidTi. The specific 
heat depends on the material details. To be more specific, 
we consider the Ni metal that has been experimentally 
investigated most extensively. 

a. Specific heat of the electrons In a free electron 
picture, the specific heat of an electron gas is Ce = 
\nn,kB{TelTF){l - 37rVlO(re/Tf)2 - ...), where is 
the electron density and Tp is the Fermi temperature^. 
However, this approximation is usually poor in the case 
of transition metals. In our model, we assume that the 
electron specific heat remains proportional to the tem- 
perature Ce = 7eT'e whcrc 7e ~ 1.5 X 10'^Jm~^K~^ is 
taken from experiments^^, which is smaller than the one 
assumed earlier—. 

h. Specific heat of lattice The phonon energy is de- 
rived from the Debyc model, Ep — J d^qhujP^nP^{T). This 
yields Ci = SNAksFoiTD/T), where Na is the Avo- 

gadro number and Fu = J^°^'^ x'^e^ /{e^ — l)'^dx is the 
Debye function. This form of the lattice specific heat is 
consistent with Pawel et ali^. 

c. Specific heat of the spins We determine the spin 
specific heat from the numerical derivative of the spin 
energy Eq. (15), as explicitly calculated in Sec. IIIII 
In Fig. [Ub) and Fig. [Ijd), we have already shown the 
temperature dependence of the spin specific heat with 
and without the external field. 



D. Summary 

We summarize below the dynamic equations that gov- 
ern the time-dependence of the three subsystem temper- 
atures after the laser pumping. 



on 



Ci{Ti) 
Cs{Ts) 



dTe_ 

dt 

<m 

dt 
dTs_ 
dt 



TeliTe, Tl) - re,(re, T,) + P(t) (29) 



Tel{T„Ti) 
^ es (^e 5 -^s) 



Tl 



(30) 
(31) 



where wc have used d£^/dt ~ Cfi(Tf^)dTf^/dt and we have 
discarded the spin-phonon interaction. In Eq. (29) we 
have inserted P{t) representing the initial laser energy 
transfer to the electrons and in Eq. (30), we have in- 
cluded a phenomenological heat diffusion of phonons to 
environment which is set at the room temperature T„„, 
this term becomes significant only at long time scale (sub- 
nanoseconds) . The functions are: 



(33) 
(34) 



(32) 



where the constants Wei and Wes are given in Eq. (24) 
and (28), and D = {6nY^^m{Ts) 



V. NUMERICAL RESULTS 

In this Section, we numerically solve our central 
Eqs. (29)-(34) for a number of plausible material pa- 
rameters. Our particular focus will be on the difference 
between our model and the previous three-temperature 
model. Since the demagnetization is mainly controlled 
by the interaction between electrons and spins, we choose 
a set of different Je^j-a large J^x representing transition 
metals (e.g., Ni, Fe and Co) and a weak J^x for some 
ferromagnetic oxides and dilute magnetic semiconduc- 
tors. Equations (29)-(31) are solved by using the follow- 
ing procedure. First, we assume that the laser instanta- 
neously heats the electron bath to Te(0) while the spin 
and lattice temperatures remain at the room tempera- 
ture Ts{0) = Ti{0) = Trm- With these initial conditions, 
we compute these temperatures after t > where the 
laser source has been turned off P{t > 0) = 0. If we only 
consider the time scale smaller than lOOps we may drop 
the heat diffusion term in Eq. (30). 

In Fig. [21 we show the typical temperature profiles after 
a low intensity laser pumping. In general, the electron- 
spin interaction is stronger than the clcctron-phonon in- 
teraction at low temperature, and spin and electron tem- 
peratures equilibrate within subpicoseconds. It takes an 
order of magnitude longer to reach the equilibrium be- 
tween lattice and the electrons. Also shown in the inset 
is the time dependent magnetization which illustrates 
the fast demagnetization and slow remagnetization. In 
Fig. [21 we show the temperature dependence of magne- 
tization for different J^x- As expected, the demagneti- 
zation time scales with the inverse of J^x while the re- 
magnetization is independent of Jex since the latter is 
controlled by the electron-lattice interaction. 

A much more interesting case is the high intensity 
of laser pumping. In this case, the spin temperature 
raises to the Curie temperature in 0.1-0.2 ps as shown 
in Fig. [31 Due to vanishingly small magnetization at the 
Curie temperature, the energy transfer between electrons 
and spins becomes negligible and thus the spin temper- 
ature stays constant for an extended period of time (a 
few ps). The electron temperature, however, continues 
to decrease due to the electron-lattice interaction which 
is not affected by the dynamic slowdown of the spins. 
Interestingly, after the electron temperature drops below 
the spin temperature, the spin system begins to heat the 
electron system and thus the electron temperature be- 
haves non-monotonically as seen in Fig. [H 

The dynamic slowdown of the spin temperature shown 
in Fig.[31is a general property of critical phenomena. Due 



Time (ps) 



Time (ps) 



FIG. 2: (Color online) Time dependence of the temperatures 
of the electrons, spins and lattice after irradiation by a low in- 
tensity laser with Te(0) = 0.7Tc, and T^O) = Tp(0) = T^™ = 
OAYTc, and Tc = 620K. The inset shows the minimum mag- 
netization (or maximum spin temperature) occurs at about 
260 femtosecond. The other parameters are: Jex = 0.15eV, 
EF = 8eV, , M/m = 10^, and ao = 0.25nm. 



FIG. 4; (Color Online) Time evolution of three temperature 
for a large laser- fluence case re(0) — 1.6. The critical slowing 
down of the spin system is identified as the plateau in the 
figure. The inset defines a slowdown time Td- The smaller 
inset shows the magnified region in the vicinity of the maxi- 
mum temperature. The other parameters are same as those 
in Fig. H 
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FIG. 3: (Color online) Time-dependent magnetization as a 
function of time in logarithmical scale for various exchange 
parameters at a fixed laser-fluence. The parameters are same 
as those of Fig. [2j 
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FIG. 5; (Color online) Log-log plot of ra versus the reduced 
temperature for Jex =0.1 and 0.2. The exponents are 5 = 
0.43 and S = 0.34 respectively. The parameters are same as 
those in Fig. [4] The dashed line is for eye-guidance. 



to the disappearance of the order parameter (the magne- 
tization m{T) in present case), the effective interaction 
reduces to zero at the critical point. In Fig. [SJ we show 
the time interval (labeled in Fig. [5]) for the critical slow- 
down as a function of the maximum spin temperature Tm 
for a given laser pumping power. As it is expected, the 
critical slowdown shows a power law, oc [1 — T„i/Tc]^^ 
with the exponent S depending on J^^ ■ In the case of very 
high intensity of the laser pumping, Tm can be very close 
to Tc and the magnetization dynamics can be extremely 
slow. In the presence of the external field, however, the 



spin system does not have a sharp phase transition any- 
more and the critical slowdown is removed, i.e., one re- 
covers the fast magnetization dynamics. In Fig. [51 we 
compare the magnetization dynamics with and without 
the magnetic field. The magnetic field suppresses the 
dynamic slowdown. 

Finally, Fig. [7] shows the exponential dependence of 
Td on the initial electron temperature which is directly 
related to the pumping laser fluence. 
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FIG. 6: (Color online) Magnetization as a function of time 
with and without the magnetic field. The parameters are 
same as those in Fig. [l] 
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FIG. 7: (Color online) Log(rd) versus the initial electron tem- 
perature for two values of Jex with (open symbols) and with- 
out (filled symbols) the magnetic field. re(0) is normalized 
initial electron temperature. 



VI. DISCUSSION 



A. Connection with experiments on LID 



tion/remagnetization time scale is governed by energy 
transfer between the baths: the demagnetization is given 
by the electron-spin interaction while remagnetization 
time is determined by the electron-phonon interactions. 

For transition metals, the electron-spin interaction is 
at least several times larger than the electron-phonon in- 
teraction. Thus, the demagnetization is faster than the 
remagnetization. For half-metals and oxidized ferromag- 
nets, the demagnetization is usually longer due to a re- 
duced electron-spin interaction. When the temperature 
increases, the demagnetization time could be significantly 
increased^; this is due to the weakening of the effective 
electron-spin interaction with a reduced magnetization. 
As the temperature approaches the Curie temperature, 
Ogasawara et alJ- observed that in all their samples, the 
demagnetization time could be enhanced by one order of 
magnitude. 

The influence of the pump intensity on the demagne- 
tization time can be similarly understood. As we have 
shown, a large pumping intensity creates high temper- 
ature electrons which heat the spin temperature to the 
Curie temperature. Thus the temperature and the pump- 
ing intensity dependence of the demagnetization involve 
the exact physics of critical slowdown. 



B. Connection with HAMR 

HAMR involves heating ferromagnets to an elevated 
temperature so that a moderate magnetic field is able to 
overcome the magnetic anisotropy for magnetization re- 
versal. Since the time scale in HAMR processes is of the 
order of nanoseconds, the dynamics studied here can be 
viewed as ultrafast, i.e., all three temperatures have al- 
ready reached equilibrium for HAMR dynamics. Even for 
the temperature close to the Curie temperature, the dy- 
namics slowdown remains "ultrafast" for HAMR as long 
as a moderate magnetic field is present. Thus, the HAMR 
dynamics could be performed in two distinct time scales: 
a fast dynamics within 10 picoseconds which determines 
the longitudinal magnetization m{T) and a slow dynam- 
ics from subnanoseconds to a few nanoseconds which de- 
termines the direction of the magnetization by the con- 
ventional Landau-Lifshitz-Gilbert equation. The detail 
calculations for HAMR dynamics will be published else- 
where. 



We now comment on the connection of our theory 
to the existing experimental results. The LID experi- 
ments performed on transition metallic fcrromagnetsii^i^ 
are usually at low laser pumping power. In these ex- 
periments, the previous phenomenological three temper- 
ature model provides an essential interpretation of de- 
magnetization: the laser induced hot electrons transfer 
their energies to spins and lattice. As discussed in Sec. 
im in our model, the demagnetization (i.e., loss of spin 
memory) occurs during the instantaneous thermalization 
of the interacting baths. Therefore, the dcmagnctiza- 



VII. CONCLUSION 

We have proposed a microscopic approach to the three 
temperature model applied to laser-induced ultrafast de- 
magnetization. The microscopic model consists of inter- 
actions among laser-excited electrons, collective spin ex- 
citations and lattice. Under the assumption of instanta- 
neous spin memory loss during the baths thermalization, 
the demagnetization problem reduces to energy transfer 
between the thcrmalizcd baths. 



10 



A self-consistent random phase approximation is de- 
veloped to model the low excitation of the spin system 
for a wide range of temperatures. A set of dynamic equa- 
tions for the time-dependent temperatures of electrons, 
spins and lattice are explicitly expressed in terms of the 
microscopic parameters. While the resulting equations 
are similar to the phcnomenological three-temperature 
model, there are important distinctions in the temper- 
ature dependent properties. In particular, the magnon 
softening plays a key role in demagnetization near Curie 
temperature where a significant slowdown of the spin dy- 
namics occurs. We have also shown that for sufHciently 
high temperatures (above the Debye temperature), the 
dynamic properties are governed by only a few parame- 
ters: the Curie and Fermi temperatures, the electron-spin 
exchange integral Jg^^, and the electron-phonon coupling 



constant Bq. The magnetization dynamic near the Curie 
temperature is rather universal. Our numerical study of 
these equations illustrates that, due to the reduction of 
the average magnetization as a function of the spin tem- 
perature, both pump intensity and sample temperature 
are responsible for a relative long demagnetization (sev- 
eral picoseconds). An external magnetic field can sup- 
press the critical dynamic slowdown. 
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